SOME IMPROVEMENTS OF ONE METHOD 
FOR PROVING INEQUALITIES BY COMPUTER 

Branko J. Malesevic 

We give some improvements of one method for proving inequalities by computer 
which is presented in the article [15]. 

Let / : [a, b] — > R be a continuous function. In this article we consider in- 
equalities in the following form: 

(1) /(*)>0 

and we give some improvements, via Propositions 1. and 3., for the method for 
proving this inequalities which is considered in the article [15]. Let us assume that 
there exist real numbers n and to such that there are finite and non-zero limits: 

(2) a= lim — and (3= lim /U '' 



x^a+ (x - a) n (b - x) m x->b- (x - a) n (b ~ x) m ' 

In the article [15] we consider only the case when n and to are non-negative integer 
points determined by: n > 1 is the multiplicity of the root x = a, otherwise n = 
if x — a is not the root; and m > 1 is the multiplicity of the root x — b, otherwise 
m = if x = b is not the root. In this case, if for the function f(x) at the point 
x = a there is an approximation of the function by Taylor polynomial of n-th order 
and at point x = b there is an approximation of the function by Taylor polynomial 
of m-th order, then [15]: 

w n\{b-a) m to! (6 -a)" 

Practically, for many examples of inequalities of the form (1), it is sufficient to 
determine limits a and /3, for the function f(x), by formulas (3) for some non- 
negative integer points n and to. Next, let us define, as it was done in [15], the 
following (^-function for the real numbers n, to and limits a, (3: 



(4) g{x) = g[ b (x) = { 



S[x) ■ xe(a,b), 



(x-a) n (b~x) m 

(3 : x = b 



The previous function is also continuous. Then the following equivalence is true: 
(5) g(x) > f{x) > 0. 
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Thus if a < or < the inequality (1) is not true. Hence, we consider only the 
cases a > and > 0. Let us notice that if the following implication is true: 

(6) g(x) > =► /(a;) > 0, 

then f(x) can have roots in the end-points of the segment [a, b]. One approach for 
proving the inequality g(x) > is based on the following statement: 

Proposition 1. Let g : [a,b] — ► R be a continuous function. Then g(x) > for 
x G [a, b] if and only if there exists a polynomial P(x) and positive real number 
S > such that: 

(7) \g(x)-P(x)\<5 
and 

(8) P(x) - S > 0. 

The necessity part of this statement is a simple consequence of the properties of 
the continuous functions and the Weierstrass approximation theorem. 

Polynomial P(x) of degree k is the minimax polynomial approximation of degree 
k of the continuous function g(x) over [a, b], if the following is true: 

(9) max \g(x) - P(x)\ < max \g(x) - P(x)\, 

for arbitrary polynomial P(x) of degree k. The minimax polynomial approximation 
P(x) of degree k is unique and the following well-known statement is true [9]: 

Theorem 2. [The Chebyshev equioscillation theorem] Polynomial P{x) is 
the minimax polynomial approximation of degree k of the continuous function g(x) 
over [a, b] if and only if there exist (fc+2) ti-points a<t < ■ ■ .<tk+i <b such that: 

(10) g(U) - P(U) = (-ITS, 5=max \g(x)-P(x)\. 

x£ [a.b] 

Polynomial Rcmez algorithm (second Rcmez algorithm [1]) is an iterative procedure 
which results are, with the chosen accuracy, the minimax polynomial approxima- 
tion and numerical estimate of the absolute error of approximation [4], [5]. The 
proof of effectiveness of polynomial Remez algorithm is given in [2] , [5] . In the case 
of two times continuously differentiable function, such that values ±8 from (10) 
exist in the end-points t$ = a, tk+i = b and in k interior ^-points also, then the rate 
of convergence of polynomial Remez algorithm is quadratic [3], [9]. On the basis 
of Maple computer program and the numapprox package, the result of polynomial 
Remez algorithm are, with the chosen accuracy, the minimax polynomial approxi- 
mation P(x) :— P(x) and numerical estimate of the absolute error 6 := 5 [11], [15]. 
If it is not possible to determine the minimax polynomial approximation in Maple 
program there appears a message that it is necessary to increase accuracy [11], [15]. 
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Next, according to the Proposition 1. for proving g(x) > it is sufficient to use 
polynomial Remez algorithm for the minimax polynomial approximation of the 
function g{x). On the other hand, it is possible to use the varieties of rational Remez 
algorithms [7], [12]. In this case we have well-known problems with convergence of 
these algorithms [10], [12]. 

On the basis of the previous considerations we can determine the method from the 
article [15] more precisely by statement: 

Proposition 3. Let f : [a, b] — > R be the continuous function for which: 

(i) there exist some real numbers n and m such that there are limits a and (3, 
which are determined by (2), as positive real numbers; 

(ii) there exist the minimax polynomial approximation P{x) and numerical es- 
timate of the absolute error 8 > for the function g{x) — gl h {x), which is 
determined by (4), such that (7) and (8) are true. 

Then f(x) > and f(x) can have roots in end-points of the segment [a, b]. 

Note 4. Let us emphasize that the previous method can be modified for other types 
of approximations of the function g{x) in order to prove inequalities g(x) > 0. 



In the article [15], as an illustration of the considered method, the following state- 
ments are proved: 

oo 

f t x — 1 

Theorem 5. For Kurepa function K(x) = / e~* — — — dt (x > 0) [8], for values 
x G [0, 1], the following inequality is true: 

(11) K(x) < K(0) x, 

where K (0) = 1.432 205 735 ... is the best possible constant. 

Remark 6. Inequality (11) is an improvement of the inequality (4.3) from Lemma 
4.2. of the article [14]. Let us emphasize that the proof of inequality (11) can be 
simply geometrically given in the following way. Let us notice that the following is 
true: 

OO OO g 

K{x)= J e -H x ^^dt>0 and K"\x)= f e-H x] ^dt>0 
o o 

for x>0. Therefore, the function K(x), for x > 0, is increasing and has only one 
inflection point c = 0.929 875 685 ... . For x e [0, c) the inequality (11) is true, 
because K(x) is concave function over [0, c). On the other hand, for x G [c, 1] the 
inequality (11) is true on the basis of: 

K(x) < K(l) = 1 < 1.331 773 289 ... = K(0) c < k'(0) x. □ 
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Theorem 7. For values x G [0, 1] the following inequality is true: 



(12) 



arc sin a; < 



tt(2 - yg) 



(VI + x -VI -a:) 



V2(4 - tt) 

7T-2V2 



+ VI + a + VI - a; 



Remark 8. Lef ws emphasize that the inequality (12) is proved in [16], on t/ie 
6asis o/ i/ie X-method Mitrinovic-Vasic [6], as an improvement of the appropriate 
inequality of L. Zhu [17]. Let us notice that the inequality (12) has been recently 
considered in [18] and [19]. 

In the article [15] the proofs of the Theorems 5. and 7. are given by the minimax 
linear approximations of the corresponding g- functions. The proofs of the Theorems 
5. and 7. are based on the assumptions that the minimax linear approximations and 
appropriate estimations of the absolute errors, given by Maple computer program, 
are correctly given under the chosen accuracy. Let us emphasize that the methods 
of the high-accuracy computations of the minimax polynomial approximations and 
the appropriate estimations of the absolute errors are considered in the Chapter 5 
of the Book [13]. Finally, let us notice that, in practice, the correctness of the 
presented method for the proving the inequality f{x) > over [a, b] is based on 
the correctness of polynomial Remez algorithm which is applied on the function 
9(x)=gl b (x) over [a, 6]. 
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